ClinPrior: an algorithm for diagnosis and novel gene discovery by network-based prioritization.

BACKGROUND: Whole-exome sequencing (WES) and whole-genome sequencing (WGS) have become indispensable tools to solve rare Mendelian genetic conditions. Nevertheless, there is still an urgent need for sensitive, fast algorithms to maximise WES/WGS diagnostic yield in rare disease patients. Most tools devoted to this aim take advantage of patient phenotype information for prioritization of genomic data, although are often limited by incomplete gene-phenotype knowledge stored in biomedical databases and a lack of proper benchmarking on real-world patient cohorts. METHODS: We developed ClinPrior, a novel method for the analysis of WES/WGS data that ranks candidate causal variants based on the patient's standardized phenotypic features (in Human Phenotype Ontology (HPO) terms). The algorithm propagates the data through an interactome network-based prioritization approach. This algorithm was thoroughly benchmarked using a synthetic patient cohort and was subsequently tested on a heterogeneous prospective, real-world series of 135 families affected by hereditary spastic paraplegia (HSP) and/or cerebellar ataxia (CA). RESULTS: ClinPrior successfully identified causative variants achieving a final positive diagnostic yield of 70% in our real-world cohort. This includes 10 novel candidate genes not previously associated with disease, 7 of which were functionally validated within this project. We used the knowledge generated by ClinPrior to create a specific interactome for HSP/CA disorders thus enabling future diagnoses as well as the discovery of novel disease genes. CONCLUSIONS: ClinPrior is an algorithm that uses standardized phenotype information and interactome data to improve clinical genomic diagnosis. It helps in identifying atypical cases and efficiently predicts novel disease-causing genes. This leads to increasing diagnostic yield, shortening of the diagnostic Odysseys and advancing our understanding of human illnesses.

RINT1 deficiency disrupts lipid metabolism and underlies a complex hereditary spastic paraplegia.

The Rad50 interacting protein 1 (Rint1) is a key player in vesicular trafficking between the ER and Golgi apparatus. Biallelic variants in RINT1 cause infantile-onset episodic acute liver failure (ALF). Here, we describe 3 individuals from 2 unrelated families with novel biallelic RINT1 loss-of-function variants who presented with early onset spastic paraplegia, ataxia, optic nerve hypoplasia, and dysmorphic features, broadening the previously described phenotype. Our functional and lipidomic analyses provided evidence that pathogenic RINT1 variants induce defective lipid-droplet biogenesis and profound lipid abnormalities in fibroblasts and plasma that impact both neutral lipid and phospholipid metabolism, including decreased triglycerides and diglycerides, phosphatidylcholine/phosphatidylserine ratios, and inhibited Lands cycle. Further, RINT1 mutations induced intracellular ROS production and reduced ATP synthesis, affecting mitochondria with membrane depolarization, aberrant cristae ultrastructure, and increased fission. Altogether, our results highlighted the pivotal role of RINT1 in lipid metabolism and mitochondria function, with a profound effect in central nervous system development.

Sphingolipid desaturase DEGS1 is essential for mitochondria-associated membrane integrity.

Sphingolipids function as membrane constituents and signaling molecules, with crucial roles in human diseases, from neurodevelopmental disorders to cancer, best exemplified in the inborn errors of sphingolipid metabolism in lysosomes. The dihydroceramide desaturase Δ4-dihydroceramide desaturase 1 (DEGS1) acts in the last step of a sector of the sphingolipid pathway, de novo ceramide biosynthesis. Defects in DEGS1 cause the recently described hypomyelinating leukodystrophy-18 (HLD18) (OMIM #618404). Here, we reveal that DEGS1 is a mitochondria-associated endoplasmic reticulum membrane-resident (MAM-resident) enzyme, refining previous reports locating DEGS1 at the endoplasmic reticulum only. Using patient fibroblasts, multiomics, and enzymatic assays, we show that DEGS1 deficiency disrupts the main core functions of the MAM: (a) mitochondrial dynamics, with a hyperfused mitochondrial network associated with decreased activation of dynamin-related protein 1; (b) cholesterol metabolism, with impaired sterol O-acyltransferase activity and decreased cholesteryl esters; (c) phospholipid metabolism, with increased phosphatidic acid and phosphatidylserine and decreased phosphatidylethanolamine; and (d) biogenesis of lipid droplets, with increased size and numbers. Moreover, we detected increased mitochondrial superoxide species production in fibroblasts and mitochondrial respiration impairment in patient muscle biopsy tissues. Our findings shed light on the pathophysiology of HLD18 and broaden our understanding of the role of sphingolipid metabolism in MAM function.

Activating cannabinoid receptor 2 preserves axonal health through GSK-3ß/NRF2 axis in adrenoleukodystrophy.

Aberrant endocannabinoid signaling accompanies several neurodegenerative disorders, including multiple sclerosis. Here, we report altered endocannabinoid signaling in X-linked adrenoleukodystrophy (X-ALD), a rare neurometabolic demyelinating syndrome caused by malfunction of the peroxisomal ABCD1 transporter, resulting in the accumulation of very long-chain fatty acids (VLCFAs). We found abnormal levels of cannabinoid receptor 2 (CB2r) and related endocannabinoid enzymes in the brain and peripheral blood mononuclear cells (PBMCs) of X-ALD patients and in the spinal cord of a murine model of X-ALD. Preclinical treatment with a selective agonist of CB2r (JWH133) halted axonal degeneration and associated locomotor deficits, along with normalization of microgliosis. Moreover, the drug improved the main metabolic disturbances underlying this model, particularly in redox and lipid homeostatic pathways, including increased lipid droplets in motor neurons, through the modulation of the GSK-3ß/NRF2 axis. JWH133 inhibited Reactive Oxygen Species elicited by excess VLCFAs in primary microglial cultures of Abcd1-null mice. Furthermore, we uncovered intertwined redox and CB2r signaling in the murine spinal cords and in patient PBMC samples obtained from a phase II clinical trial with antioxidants (NCT01495260). These findings highlight CB2r signaling as a potential therapeutic target for X-ALD and perhaps other neurodegenerative disorders that present with dysregulated redox and lipid homeostasis.

Epigenetic profiling linked to multisystem inflammatory syndrome in children (MIS-C): A multicenter, retrospective study.

Background: Most children and adolescents infected with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remain asymptomatic or develop a mild coronavirus disease 2019 (COVID-19) that usually does not require medical intervention. However, a small proportion of pediatric patients develop a severe clinical condition, multisystem inflammatory syndrome in children (MIS-C). The involvement of epigenetics in the control of the immune response and viral activity prompted us to carry out an epigenomic study to uncover target loci regulated by DNA methylation that could be altered upon the appearance of MIS-C. Methods: Peripheral blood samples were recruited from 43 confirmed MIS-C patients. 69 non-COVID-19 pediatric samples and 15 COVID-19 pediatric samples without MIS-C were used as controls. The cases in the two groups were mixed and divided into discovery (MIS-C = 29 and non-MIS-C = 56) and validation (MIS-C = 14 and non-MIS-C = 28) cohorts, and balanced for age, gender and ethnic background. We interrogated 850,000 CpG sites of the human genome for DNA methylation variants. Findings: The DNA methylation content of 33 CpG loci was linked with the presence of MIS-C. Of these sites, 18 (54.5%) were located in described genes. The top candidate gene was the immune T-cell mediator ZEB2; and others highly ranked candidates included the regulator of natural killer cell functional competence SH2D1B; VWA8, which contains a domain of the Von Willebrand factor A involved in the pediatric hemostasis disease; and human leukocyte antigen complex member HLA-DRB1; in addition to pro-inflammatory genes such as CUL2 and AIM2. The identified loci were used to construct a DNA methylation profile (EPIMISC) that was associated with MIS-C in both cohorts. The EPIMISC signature was also overrepresented in Kawasaki disease patients, a childhood pathology with a possible viral trigger, that shares many of the clinical features of MIS-C. Interpretation: We have characterized DNA methylation loci that are associated with MIS-C diagnosis. The identified genes are likely contributors to the characteristic exaggerated host inflammatory response observed in these patients. The described epigenetic signature could also provide new targets for more specific therapies for the disorder. Funding: Unstoppable campaign of Josep Carreras Leukaemia Foundation, Fundació La Marató de TV3, Cellex Foundation and CERCA Programme/Generalitat de Catalunya.

Accelerated biological aging in COVID-19 patients.

Chronological age is a risk factor for SARS-CoV-2 infection and severe COVID-19. Previous findings indicate that epigenetic age could be altered in viral infection. However, the epigenetic aging in COVID-19 has not been well studied. In this study, DNA methylation of the blood samples from 232 healthy individuals and 413 COVID-19 patients is profiled using EPIC methylation array. Epigenetic ages of each individual are determined by applying epigenetic clocks and telomere length estimator to the methylation profile of the individual. Epigenetic age acceleration is calculated and compared between groups. We observe strong correlations between the epigenetic clocks and individual's chronological age (r > 0.8, p < 0.0001). We also find the increasing acceleration of epigenetic aging and telomere attrition in the sequential blood samples from healthy individuals and infected patients developing non-severe and severe COVID-19. In addition, the longitudinal DNA methylation profiling analysis find that the accumulation of epigenetic aging from COVID-19 syndrome could be partly reversed at late clinic phases in some patients. In conclusion, accelerated epigenetic aging is associated with the risk of SARS-CoV-2 infection and developing severe COVID-19. In addition, the accumulation of epigenetic aging from COVID-19 may contribute to the post-COVID-19 syndrome among survivors.

Biallelic variants in SLC35B2 cause a novel chondrodysplasia with hypomyelinating leukodystrophy.

Sulphated proteoglycans are essential in skeletal and brain development. Recently, pathogenic variants in genes encoding proteins involved in the proteoglycan biosynthesis have been identified in a range of chondrodysplasia associated with intellectual disability. Nevertheless, several patients remain with unidentified molecular basis. This study aimed to contribute to the deciphering of new molecular bases in patients with chondrodysplasia and neurodevelopmental disease. Exome sequencing was performed to identify pathogenic variants in patients presenting with chondrodysplasia and intellectual disability. The pathogenic effects of the potentially causative variants were analysed by functional studies. We identified homozygous variants (c.1218_1220del and c.1224_1225del) in SLC35B2 in two patients with pre- and postnatal growth retardation, scoliosis, severe motor and intellectual disabilities and hypomyelinating leukodystrophy. By functional analyses, we showed that the variants affect SLC35B2 mRNA expression and protein subcellular localization leading to a functional impairment of the protein. Consistent with those results, we detected proteoglycan sulphation impairment in SLC35B2 patient fibroblasts and serum. Our data support that SLC35B2 functional impairment causes a novel syndromic chondrodysplasia with hypomyelinating leukodystrophy, most likely through a proteoglycan sulphation defect. This is the first time that SLC35B2 variants are associated with bone and brain development in human.

Diagnosis of Genetic White Matter Disorders by Singleton Whole-Exome and Genome Sequencing Using Interactome-Driven Prioritization.

BACKGROUND AND OBJECTIVES: Genetic white matter disorders (GWMD) are of heterogeneous origin, with >100 causal genes identified to date. Classic targeted approaches achieve a molecular diagnosis in only half of all patients. We aimed to determine the clinical utility of singleton whole-exome sequencing and whole-genome sequencing (sWES-WGS) interpreted with a phenotype- and interactome-driven prioritization algorithm to diagnose GWMD while identifying novel phenotypes and candidate genes. METHODS: A case series of patients of all ages with undiagnosed GWMD despite extensive standard-of-care paraclinical studies were recruited between April 2017 and December 2019 in a collaborative study at the Bellvitge Biomedical Research Institute (IDIBELL) and neurology units of tertiary Spanish hospitals. We ran sWES and WGS and applied our interactome-prioritization algorithm based on the network expansion of a seed group of GWMD-related genes derived from the Human Phenotype Ontology terms of each patient. RESULTS: We evaluated 126 patients (101 children and 25 adults) with ages ranging from 1 month to 74 years. We obtained a first molecular diagnosis by singleton WES in 59% of cases, which increased to 68% after annual reanalysis, and reached 72% after WGS was performed in 16 of the remaining negative cases. We identified variants in 57 different genes among 91 diagnosed cases, with the most frequent being RNASEH2B, EIF2B5, POLR3A, and PLP1, and a dual diagnosis underlying complex phenotypes in 6 families, underscoring the importance of genomic analysis to solve these cases. We discovered 9 candidate genes causing novel diseases and propose additional putative novel candidate genes for yet-to-be discovered GWMD. DISCUSSION: Our strategy enables a high diagnostic yield and is a good alternative to trio WES/WGS for GWMD. It shortens the time to diagnosis compared to the classical targeted approach, thus optimizing appropriate management. Furthermore, the interactome-driven prioritization pipeline enables the discovery of novel disease-causing genes and phenotypes, and predicts novel putative candidate genes, shedding light on etiopathogenic mechanisms that are pivotal for myelin generation and maintenance.

Biallelic PI4KA variants cause a novel neurodevelopmental syndrome with hypomyelinating leukodystrophy.

Phosphoinositides are lipids that play a critical role in processes such as cellular signalling, ion channel activity and membrane trafficking. When mutated, several genes that encode proteins that participate in the metabolism of these lipids give rise to neurological or developmental phenotypes. PI4KA is a phosphoinositide kinase that is highly expressed in the brain and is essential for life. Here we used whole exome or genome sequencing to identify 10 unrelated patients harbouring biallelic variants in PI4KA that caused a spectrum of conditions ranging from severe global neurodevelopmental delay with hypomyelination and developmental brain abnormalities to pure spastic paraplegia. Some patients presented immunological deficits or genito-urinary abnormalities. Functional analyses by western blotting and immunofluorescence showed decreased PI4KA levels in the patients' fibroblasts. Immunofluorescence and targeted lipidomics indicated that PI4KA activity was diminished in fibroblasts and peripheral blood mononuclear cells. In conclusion, we report a novel severe metabolic disorder caused by PI4KA malfunction, highlighting the importance of phosphoinositide signalling in human brain development and the myelin sheath.

Epigenome-wide association study of COVID-19 severity with respiratory failure.

BACKGROUND: Patients infected with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for the coronavirus disease 2019 (COVID-19), exhibit a wide spectrum of disease behaviour. Since DNA methylation has been implicated in the regulation of viral infections and the immune system, we performed an epigenome-wide association study (EWAS) to identify candidate loci regulated by this epigenetic mark that could be involved in the onset of COVID-19 in patients without comorbidities. METHODS: Peripheral blood samples were obtained from 407 confirmed COVID-19 patients ≤ 61 years of age and without comorbidities, 194 (47.7%) of whom had mild symptomatology that did not involve hospitalization and 213 (52.3%) had a severe clinical course that required respiratory support. The set of cases was divided into discovery (n = 207) and validation (n = 200) cohorts, balanced for age and sex of individuals. We analysed the DNA methylation status of 850,000 CpG sites in these patients. FINDINGS: The DNA methylation status of 44 CpG sites was associated with the clinical severity of COVID-19. Of these loci, 23 (52.3%) were located in 20 annotated coding genes. These genes, such as the inflammasome component Absent in Melanoma 2 (AIM2) and the Major Histocompatibility Complex, class I C (HLA-C) candidates, were mainly involved in the response of interferon to viral infection. We used the EWAS-identified sites to establish a DNA methylation signature (EPICOVID) that is associated with the severity of the disease. INTERPRETATION: We identified DNA methylation sites as epigenetic susceptibility loci for respiratory failure in COVID-19 patients. These candidate biomarkers, combined with other clinical, cellular and genetic factors, could be useful in the clinical stratification and management of patients infected with the SARS-CoV-2. FUNDING: The Unstoppable campaign of the Josep Carreras Leukaemia Foundation, the Cellex Foundation and the CERCA Programme/Generalitat de Catalunya.

Neutralizing Autoantibodies to Type I IFNs in >10% of Patients with Severe COVID-19 Pneumonia Hospitalized in Madrid, Spain.

BACKGROUND: In a recent study, autoantibodies neutralizing type I interferons (IFNs) were present in at least 10% of cases of critical COVID-19 pneumonia. These autoantibodies neutralized most type I IFNs but rarely IFN-beta. OBJECTIVES: We aimed to define the prevalence of autoantibodies neutralizing type I IFN in a cohort of patients with severe COVID-19 pneumonia treated with IFN-beta-1b during hospitalization and to analyze their impact on various clinical variables and outcomes. METHODS: We analyzed stored serum/plasma samples and clinical data of COVID-19 patients treated subcutaneously with IFN-beta-1b from March to May 2020, at the Infanta Leonor University Hospital in Madrid, Spain. RESULTS: The cohort comprised 47 COVID-19 patients with severe pneumonia, 16 of whom (34%) had a critical progression requiring ICU admission. The median age was 71 years, with 28 men (58.6%). Type I IFN-alpha- and omega-neutralizing autoantibodies were found in 5 of 47 patients with severe pneumonia or critical disease (10.6%), while they were not found in any of the 118 asymptomatic controls (p = 0.0016). The autoantibodies did not neutralize IFN-beta. No demographic, comorbidity, or clinical differences were seen between individuals with or without autoantibodies. We found a significant correlation between the presence of neutralizing autoantibodies and higher C-reactive protein levels (p = 5.10e-03) and lower lymphocyte counts (p = 1.80e-02). No significant association with response to IFN-beta-1b therapy (p = 0.34) was found. Survival analysis suggested that neutralizing autoantibodies may increase the risk of death (4/5, 80% vs 12/42, 28.5%). CONCLUSION: Autoantibodies neutralizing type I IFN underlie severe/critical COVID-19 stages in at least 10% of cases, correlate with increased C-RP and lower lymphocyte counts, and confer a trend towards increased risk of death. Subcutaneous IFN-beta treatment of hospitalized patients did not seem to improve clinical outcome. Studies of earlier, ambulatory IFN-beta treatment are warranted.

Shorter androgen receptor polyQ alleles protect against life-threatening COVID-19 disease in European males.

BACKGROUND: While SARS-CoV-2 similarly infects men and women, COVID-19 outcome is less favorable in men. Variability in COVID-19 severity may be explained by differences in the host genome. METHODS: We compared poly-amino acids variability from WES data in severely affected COVID-19 patients versus SARS-CoV-2 PCR-positive oligo-asymptomatic subjects. FINDINGS: Shorter polyQ alleles (≤22) in the androgen receptor (AR) conferred protection against severe outcome in COVID-19 in the first tested cohort (both males and females) of 638 Italian subjects. The association between long polyQ alleles (≥23) and severe clinical outcome (p = 0.024) was also validated in an independent cohort of Spanish men <60 years of age (p = 0.014). Testosterone was higher in subjects with AR long-polyQ, possibly indicating receptor resistance (p = 0.042 Mann-Whitney U test). Inappropriately low serum testosterone level among carriers of the long-polyQ alleles (p = 0.0004 Mann-Whitney U test) predicted the need for intensive care in COVID-19 infected men. In agreement with the known anti-inflammatory action of testosterone, patients with long-polyQ and age ≥60 years had increased levels of CRP (p = 0.018, not accounting for multiple testing). INTERPRETATION: We identify the first genetic polymorphism that appears to predispose some men to develop more severe disease. Failure of the endocrine feedback to overcome AR signaling defects by increasing testosterone levels during the infection leads to the polyQ tract becoming dominant to serum testosterone levels for the clinical outcome. These results may contribute to designing reliable clinical and public health measures and provide a rationale to test testosterone as adjuvant therapy in men with COVID-19 expressing long AR polyQ repeats. FUNDING: MIUR project "Dipartimenti di Eccellenza 2018-2020" to Department of Medical Biotechnologies University of Siena, Italy (Italian D.L. n.18 March 17, 2020) and "Bando Ricerca COVID-19 Toscana" project to Azienda Ospedaliero-Universitaria Senese. Private donors for COVID-19 research and charity funds from Intesa San Paolo.

Case Report: Benign Infantile Seizures Temporally Associated With COVID-19.

Background: Non-febrile illness seizures may present in previously healthy children as afebrile seizures associated with minor infections, such as mild gastroenteritis or respiratory tract infections, and are linked to a genetic predisposition. For the novel human coronavirus SARS-CoV-2, causing COVID-19, fever, cough, and gastrointestinal complaints are the most common symptoms in children, and a hyperimmune response may be present. No detailed temporally associated neurological complications have been documented in pediatric case series so far. Case description: We present the case of a 3-months-old girl with non-febrile repeated seizures in a COVID-19 family setting. The infant started with a mild fever and cough that lasted for 2 days. At day 6 from onset, the girl presented with two focal motor seizures with impaired consciousness and awareness. All investigations ruled out signs of meningo-encephalitis or active epilepsy, including normal electroencephalogram and cerebral magnetic resonance imaging. PCR from nasal and throat swabs was positive for SARS-CoV-2. Remarkably, blood ferritin and D-dimer levels were increased. At day 9, the infant presented another afebrile motor seizure, and levetiracetam dose was modified there was a favorable response within 3 months of the follow-up. Much interest has been raised with regards to host genetic determinants to disease severity and susceptibility to COVID-19. We thus performed whole exome sequencing, revealing a pathogenic frameshift mutation in the PRRT2 gene in both the mother and the infant. The mother had presented two late infantile febrile convulsions with normal outcome afterwards. Discussion: The hyperimmune response described in adult cases with COVID-19 can be seen in infants, even in the absence of respiratory symptoms. Moreover, COVID-19 may present in infants as non-febrile seizures, triggering early onset seizures in infants with a genetic predisposition. In this pandemic situation, precision medicine using massive sequencing can shed light on underlying molecular mechanisms driving the host response to COVID-19.

HNRNPH1-related syndromic intellectual disability: Seven additional cases suggestive of a distinct syndromic neurodevelopmental syndrome.

Pathogenic variants in HNRNPH1 were first reported in 2018. The reported individual, a 13 year old boy with a c.616C>T (p.R206W) variant in the HNRNPH1 gene, was noted to have overlapping symptoms with those observed in HNRNPH2-related X-linked intellectual disability, Bain type (MRXSB), specifically intellectual disability and dysmorphic features. While HNRNPH1 variants were initially proposed to represent an autosomal cause of MRXSB, we report an additional seven cases which identify phenotypic differences from MRXSB. Patients with HNRNPH1 pathogenic variants diagnosed via WES were identified using clinical networks and GeneMatcher. Features unique to individuals with HNRNPH1 variants include distinctive dysmorphic facial features; an increased incidence of congenital anomalies including cranial and brain abnormalities, genitourinary malformations, and palate abnormalities; increased incidence of ophthalmologic abnormalities; and a decreased incidence of epilepsy and cardiac defects compared to those with MRXSB. This suggests that pathogenic variants in HNRNPH1 result in a related, but distinct syndromic cause of intellectual disability from MRXSB, which we refer to as HNRNPH1-related syndromic intellectual disability.

A deep intronic splice variant advises reexamination of presumably dominant SPG7 Cases.

OBJECTIVE: To identify causative mutations in a patient affected by ataxia and spastic paraplegia. METHODS: Whole-exome sequencing (WES) and whole-genome sequencing (WGS) were performed using patient's DNA sample. RT-PCR and cDNA Sanger sequencing were performed on RNA extracted from patient's fibroblasts, as well as western blot. RESULTS: A novel missense variant in SPG7 (c.2195T> C; p.Leu732Pro) was first found by whole-exome sequencing (WES), while the second, also unreported, deep intronic variant (c.286 + 853A>G) was identified by whole-genome sequencing (WGS). RT-PCR confirmed the in silico predictions showing that this variant activated a cryptic splice site, inducing the inclusion of a pseudoexon into the mRNA sequence, which encoded a premature stop codon. Western blot showed decreased SPG7 levels in patient's fibroblasts. INTERPRETATION: Identification of a deep intronic variant in SPG7, which could only have been detected by performing WGS, led to a diagnosis in this HSP patient. This case challenges the notion of an autosomal dominant inheritance for SPG7, and illustrates the importance of performing WGS subsequently or alternatively to WES to find additional mutations, especially in patients carrying one variant in a gene causing a predominantly autosomal recessive disease.